SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 

PHYSICS 

STATISTICAL PHYSICS, 
THERMODYNAMICS AND 
NONLINEAR DYNAMICS 

PHYSICAL SCIENCES 

COMPLEX NETWORKS 



Self-organization versus top-down 
planning in the evolution of a city 

Marc Barthelemy 12 , Patricia Bordin 3 4 , Henri Berestycki 2 &Maurizio Gribaudi 5 

1 1nstitut de Physique Theorique, CEA, CNRS-URA 2306, F-91 1 91 , Gif-sur-Yvette, France, 2 Centre d'Analyse et de Mathematiques 
Sociales, EHESS, 1 90-1 98 avenue de France, 75244 Paris, 3 Universite Paris-Est, Institut de Recherche en Constructibilite, ESTP, F- 
94230, Cachan, France, 4 Universite Paris Diderot, Sorbonne Paris Cite, Institut des Energies de Demain (IED), 75205 Paris, France, 
5 Laboratoire de Demographie et Histoire Sociale, EHESS, 1 90-1 98 avenue de France, 75244 Paris, France. 



Received 
26 April 201 3 

Accepted 
19 June 2013 

Published 
8 July 2013 



Correspondence and 
requests for materials 
should be addressed to 
M.B. (marc. 
barthelemy@cea.fr) 



Interventions of central, top-down planning are serious limitations to the possibility of modelling the 
dynamics of cities. An example is the city of Paris (France), which during the 19th century experienced large 
modifications supervised by a central authority, the 'Haussmann period'. In this article, we report an 
empirical analysis of more than 200 years (1789-2010) of the evolution of the street network of Paris. We 
show that the usual network measures display a smooth behavior and that the most important quantitative 
signatures of central planning is the spatial reorganization of centrality and the modification of the block 
shape distribution. Such effects can only be obtained by structural modifications at a large-scale level, with 
the creation of new roads not constrained by the existing geometry. The evolution of a city thus seems to 
result from the superimposition of continuous, local growth processes and punctual changes operating at 
large spatial scales. 

A city is a highly complex system where a large number of agents interact, leading to a dynamics seemingly 
difficult to understand. Many studies in history, geography, spatial economics, sociology, or physics 
discuss various facets of the evolution of the city 1 " 10 . From a very general perspective, the large number 
and the diversity of agents operating simultaneously in a city suggest the intriguing possibility that cities are an 
emergent phenomenon ruled by self- organization 2 . On the other hand, the existence of central planning inter- 
ventions might minimize the importance of self- organization in the course of evolution of cities. Central planning 
-here understood as a top-down process controlled by a central authority - plays an important role in the city, 
leaving long standing traces, even if the time horizon of planners is limited and much smaller than the age of the 
city. One is thus confronted with the question of the possiblity of modelling a city and its expansion as a self- 
organized phenomenon. Indeed central planning could be thought of as an external perturbation, as if it were 
foreign to the self- organized development of a city. The recent digitization and georeferentiation of old maps will 
enable us to test quantitatively this effect, at least at the level of the structure of the road network. Such a 
transportation network is a crucial ingredient in cities as it allows individuals to work, transport and exchange 
goods, etc., and the evolution of this network reflects the evolution of the population and activity densities 1112 . 
These network aspects were first studied in the 1960s in quantitative geography 13 , and in the last decade, complex 
networks theory has provided significant contributions to the quantitative characterization of urban street 
patterns 14 " 24 . 

In this article, we will consider the case of the evolution of the street network of Paris over more than 200 years 
with a particular focus on the 19th century, period when Paris experienced large transformations under the 
guidance of Baron Haussmann 25 . It would be difficult to describe the social, political, and urbanistic importance 
and impact of Haussmann works in a few lines here and we refer the interested reader to the existing abundant 
literature on the subject (see 26 , and 25 and references therein). Essentially, until the middle of the 19th century, 
central Paris has a medieval structure composed of many small and crowded streets, creating congestion and, 
according to some contemporaries, probably health problems. In 1852, Napoleon III commissioned Haussmann 
to modernize Paris by building safer streets, large avenues connected to the new train stations, central or symbolic 
squares (such as the famous place de l'Etoile, place de la Nation and place du Pantheon), improving the traffic flow 
and, last but not least, the circulation of army troops. Haussmann also built modern housing with uniform 
building heights, new water supply and sewer systems, new bridges, etc (see Fig. 1 where we show how dramatic 
the impact of Haussmann transformations are). The case of Paris under Haussmann provides an interesting 
example where changes due to central planning are very important and where a naive modelling is bound to fail. 
We analyze here in detail the effect of these planned transformations on the street network. By introducing 
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Figure 1 | Illustration on a small area of the impact of Haussmann's 
transformations. On the yellow background, we show the parcel 
distribution before Haussmann (extracted from the Vasserot cadastre, 
1808-1836), and in brown we show the new buildings delineating the new 
streets as designed by Haussmann and as they appeared in 1 888. We can see 
on this example that the Haussmann plan implied a large number of 
destruction and rebuilding: approximately 28, 000 houses were destroyed 
and 100, 000 were built 25 (figure created from authors' data with 
QuantumGIS software: www.qgis.org). 

physical quantitative measures associated with this network, we are 
able to compare the effect of the Hausmann transformation of the 
city with its 'natural' evolution over other periods. 

By digitizing historical maps (for details on the sources used to 
construct the maps, see the Methods section) into a Geographical 
Information System (GIS) environment, we reconstruct the detailed 
road system (including minor streets) at six different moments in 
time, t = 1, 2, 6, respectively corresponding to years: 1789, 1826, 
1836, 1888, 1999, 2010. For each time, we constructed the associated 
primal graph G t (see the Methods section and 23 ' 24 ), i.e. the graph 
where the nodes represent street junctions and the links correspond 
to road segments. In particular, it is important to note that we have 
thus snapshots of the street network before Haussmann works 
(1789-1836) and after (1888-2010). This allows us to study quanti- 
tatively the effect of such central planning. 

In Fig. 2(a), we display the map of Paris as it was in 1789 on top of 
the current map (2010). In order to use a single basis for comparison, 
we limited our study over time to the portion corresponding to 1789. 
We note here that the evolution of the outskirts and small villages in 
the surroundings has certainly an impact on the evolution of Paris 
and even if we focus here (mainly because of data availability reasons) 



on the structural modifications of the inner structure of Paris, a study 
at a larger scale will certainly be needed for capturing the whole 
picture of the evolution of this city. We then have 6 maps for different 
times and for the same area (of order 34 km 2 ). We also represent on 
Fig. 2(b), the new streets created during the Haussmann period 
which covers roughly the second half of the 19th century. Even if 
we observe some evolution outside of this portion, most of the 
Haussmann works are comprised within this portion. 

Results 

Simple measures. In the following we will study the structure of the 
graph G t at different times t (see the Methods section for precise 
definitions), having in mind that our goal is to identify the most 
important quantitative signatures of central planning during the 
evolution of this road network. 

First basic measures include the evolution of the number of nodes 
N, edges £, and total length L tot of the networks (restricted to the area 
corresponding to 1789). In Fig. 3 we show the results for these indi- 
cators which display a clear acceleration during the Haussmann 
period (1836-1888). The number of nodes increased from about 
3000 in 1836 to about 6000 in 1888 and the total length increase 
from about 400 kms to almost 700 kms, all this in about 50 years. It is 
interesting to note that this node increase corresponds essentially to 
an important increase in the population. In particular, we note (see 
the Supplementary Information for more details) that the number of 
nodes AT is proportional to the population P and that the correspond- 
ing increase rate is of order dN/dP ~ 0.0021, similar to what was 
measured in a previous study about a completely different area 24 . The 
rapid increase of nodes during the Haussmann period is thus largely 
due to demographic pressure. Now, if we want to exclude exogeneous 
effects and focus on the structure of networks, we can plot the various 
indicators such as the number of edges and the total length versus the 
number of nodes taken as a time clock. The results shown Fig. 3(d-f) 
display a smoother behavior. In particular, £ is a linear function of N, 
demonstrating that the average degree is essentially constant (k) ~ 
3.0 since 1789. The total length versus N also displays a smooth 
behavior consistent with a perturbed lattice 23 . Indeed, if the segment 
length € x is roughly constant and equal to t\ = \j yj~p where p = N/A 
is the density of nodes (A is the area considered here), we then obtain 
for the total length 

l w = { ^-Van (i) 

A fit of the type ay^N is shown in Fig. 3(d) and the value of a 
measured gives an estimate of the area A ~ 29. 7km 2 , in agreement 
with the actual value A = 33.6 km 2 (for the 1789 portion). This 
agreement demonstrates that all the networks at different times are 
not far from a perturbed lattice. 




Figure 2 | (a) Map of Paris in 1789 superimposed on the map of current 2010 Paris. In the whole study, we focus on the Haussmann modifications 
and limited ourselves to the 1789 portion of the street network, (b) Map of Haussmann modifications. The grey lines represent the road network 
in 1836, the green lines represent the Haussmann modifications which are basically all contained in the 1789 area (figure created from authors' data with 
ArcGiS 10.0 (ESRI): www.esri.com/software/arcgis/). 
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Figure 3 | Top panels: Number of (a) nodes, (b) total length (kms), and (c) rescaled average route distance versus time. Bottom panels: Number of 
(d) edges, (e) total length (kms), and (f) the rescaled average route distance versus the number of nodes N. In (d) the dashed (blue) line is a linear 
fit with slope 1.55 (r 2 = 0.99) consistent with constant average degree of order (A:) ~ 3, and in (e) the dashed (green) line a square root fit of the form ay/N 
withfl= 8.44 kmsfr 2 = 0.99). Based on a perturbed lattice picture this gives an area equal to A ~ 29.7km 2 consistent with the actual value (A= 33.6 km 2 ). 
In (f), we show the rescaled average shortest route versus AT which decreases showing that the denser the network and the easier it is to navigate 
from one node to the other (if delays at junctions are neglected). 



We also plot the average route distance d R defined as the average 
over all pairs of nodes of the shortest route between them (see 
Methods for more details). For a two dimensional spatial network, 
we expect this quantity to scale as cIr ~ y/N and thus increases with 
N The ratio d R j y/N is thus better suited to measure the efficiency of 
the network and we observe (Fig. 3(c,f)) that it decreases with time 
and N. This result simply demonstrate that if we neglect delays at 
junctions, it becomes easier to navigate in the network as it gets 
denser. 

Typology of new links. We can have three different types of new 
links depending on the number of new nodes they connect. We 
denote by E { (i = 1, 2, 3) the number of new links appearing at 
time t + 1 connecting i new nodes. For example E 0 counts the new 
links appearing at time t + 1 connecting two nodes existing at time t. 
In order to categorize more precisely these new links, we use the 
betweenness centrality impact d defined in 24 and which measures 
how a new link (absent at time t and present at time t + 1) affects 
the average betweenness centrality (see Methods section for 
definitions of the betweenness centrality impact 3). In 24 , the 
distribution of this quantity displays two peaks which corresponds 
to two types of links belonging to two distinct processes: densification 
and exploration 24 . We first observe (see Figure 2 of SI) that in the first 
period, the majority of new links are of the E 2 type and correspond to 
construction of new streets with new nodes. We see that the 
Haussmann transition period (1836-1888) is not particularly 
different from the other previous periods. In the modern period 
(after 1999), E 0 becomes dominant and consistent with the idea of 
a mature street network where densification dominates the evolution 
of the urban tissue. Obviously, this is also an effect of limiting 
ourselves to the 1789 portion: in a wider area, many new roads 
were created and both densification and exploration coexist. We 
note here that the structure of the street network of central Paris 
remained remarkably stable from 1888 until now (and in this 
period also, densification was the main process in this area). 



We then plot the distribution of this quantity 3 for the different 
transition periods and the result is shown in Fig. 4. These figures 
show that for all periods most new links belong to the densification 
process with a small peak of exploration in the period 1836-1888. 
In well -developed, mature systems, it is expected that densification 
is the dominant growth mechanism. Here also, we see that the 
Haussmann period is not significantly different from previous 
periods. 

Evolution of the spatial distribution of centrality. The betweenness 
centrality (BC) g v (i) of a node i is defined in the Methods section and 
essentially measures the fraction of times a given node is used in the 
shortest paths connecting any pair of nodes in the network, and is 
thus a measure of the contribution of a link in the organisation of 
flows in the network 27 . In our case where we consider a limited 
portion of a spatial network, two important effects need to be 
taken into consideration. First, as we consider a portion, only 
paths within this portion are taken into account in the calculation 
of the BC and this usually does not reflect the reality of the actual 
origin -destination matrix. In particular, flows with the exterior of the 
portion and surrounding villages are not taken into account. As a 
result, the BC will be able to detect important routes and nodes in 
the internal structure of the network but will miss large-scale 
communication roads such as a north-south or a East- West road 
connecting the portion with the surroundings of Paris. In 24 , the 
scale of the network was large enough so that the BC could recover 
important central roads such as Roman streets. The BC in the present 
case has then to be used as a structural probe of the network, enabling 
us to track the important modifications. The second point concerns 
the spatial distribution of the BC which will be important in the 
following. For a lattice the most central nodes (see the discussion 
in 23 for example) are close to the barycenter of the nodes: spatial 
centrality and betweenness centrality are then usually strongly 
correlated. In 16 and 17 it is shown that the most central points 
display interesting spatial structures which still need to be 
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Figure 4 | Betweenness centrality impact distribution for the periods 1789 -> 1826, 1826 -> 1836, 1836 -> 1888, 1888 -> 1999, 1999 -> 2010. 

This figure shows that densification is the main process for this portion of Paris and that from this point of view, the Haussmann period seems to be rather 
smooth and comparable to other periods. 



understood, but which represent an important signature of the 
networks' topology. 

We first consider the time evolution of the node betweenness 
centrality (with similar results for the edge BC). In the SI (see figure 3 
of SI), we show the distribution of the node BC at different times. 
Apart from the fact that the average BC varies, we see that the tail of 
the distribution remains constant in time, showing that the statistics 
of very central nodes is not modified. From this point of view, the 
evolution of the road network follows a smooth behavior, even in the 
Haussmann period. 

So far, most of the measures indicate that the evolution of the street 
network follows simple densification and exploration rules and is 
very similar to other areas studied 24 . At this point, it appears that 
Haussmann works didn't change radically the structure of the city. 
However, we can suspect that Haussmann's impact is very important 
on congestion and traffic and should therefore be seen on the spatial 



distribution of centrality. In the figure 5, we show the maps of Paris at 
different times and we indicate the most central nodes (such that 
their centrality g v (i) is larger than maxg v /a with a = 10 see the SI, for 
a discussion on the effect of the value of a). We can clearly see here 
that the spatial distribution of the BC is not stable, displays large 
variations, and is not uniformly distributed over the Paris area (we 
represented here the node centrality, and similar results are obtained 
for the edge centrality, see the SI for plots for the edge centrality and 
more details). In particular, we see that between 1836 and 1888, the 
Haussmann works had a dramatical impact on the spatial structure 
of the centrality, especially near the heart of Paris. Central roads 
usually persist in time 24,28 , but in our case, the Haussmann reorgan- 
ization was acting precisely at this level by redistributing the shortest 
paths which had certainly an impact on congestion inside the city. 
After Haussmann we observe a large stability of the network until 
nowadays. 




Figure 5 | Spatial distribution of the most central nodes (with centrality g v such that g v > max g v /10). We observe for the different periods important 
reorganizations of the spatial distribution of centrality, corresponding to different specific interventions. In particular, we observe a very important 
redistribution of centrality during the Haussmann period with the appearance of a reticulated structure on the 1888 map. 
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Figure 6 | Top panels: (a) number of nodes in a disk of radius r from the barycenter of Paris and (b) its variation versus r. As expected the largest 
variation occurred during the Haussmann period. Bottom panels: (c) Number of nodes at distance r and with centrality larger than go (go = max gv/10) 
and (d) its variation. The thick green line in the right panels indicate the Haussmann transition 1836-1888. We see here that during the Haussmann 
period (and also in the 1789-1826 period), there is a large decrease of the number of central nodes in the central region of Paris (r < 2,000 meters). 



It is interesting to note that these maps also provide details about 
the evolution of the road network of Paris during other periods which 
seems to reflect what happened in reality and which we can relate to 
specific local interventions. For example, in the period 1789-1826 
between the French Revolution and the Napoleonic empire, the maps 
shown in Fig. 5 display large variations with redistribution of central 
nodes which probably reflects the fact that many religious and aris- 
tocratic domains and properties were sold and divided in order to 
create new houses and new roads, improving congestion inside Paris. 
During the period 1826-1836 which corresponds roughly to the 
beginning of the the July Monarchy, the maps in Fig. 5 suggests an 
important reorganization on the east side of Paris. This seems to 
correspond very well to the creation during that period of a new 
channel in this area (the channel 'Saint Martin') which triggered 
many transformations in the eastern part of the network. 

In order to analyse the spatial redistribution effect more quantita- 
tively, we compute various quantities inside a disk of radius r centered 
on the barycenter of all nodes (which stays approximately at the same 
location in time). We first study the number of nodes N(r) (Fig. 6), its 
variation 3N(r) between t and t + 1, and the number of central nodes 
(such that g v (t) > max g v /!0). We see that the largest variation of the 
number of nodes (see Fig. 6b) is indeed in the Haussmann period 
1836-1888, especially for distance r > 1,500 meters. More interesting, 
is the variation of the most central nodes (Fig. 6d). In particular, we 
observe that during the pre- Haussmann period, even if in the period 
1789-1826 there was an improvement of centrality concentration, 
there is an accumulation of central nodes both at short distances (r 
< 2,500 meters) and at long distances (r > 2,500 meters) in the 
following period (1826-1836). As a result, visually clear in Fig. 5, 
there is a large concentration of centrality in the center of Paris until 
1836 at least. The natural consequence of this concentration is that 
the center of Paris was very probably very congested at that time. In 
this respect, what happens under the Haussmann supervision is nat- 
ural as he acts on the spatial organization of centrality. We see indeed 
that in 1888, the most central nodes form a more reticulated structure 
excluding concentration of centrality. A structure which remained 
stable until now. Interestingly, we note that Haussmann's new roads 
and avenues represent approximately 6% of the total length only 



(compared to nowadays network), which is a small fraction, consid- 
ered that it has a very important impact on the centrality spatial 
organization. 

This reorganization of centrality was undertaken with creation of 
new roads and avenues destroying parts of the original pattern (see 
Fig. 1 and Fig. 2(b)) resulting in the modification of the geometrical 
structure of blocks (defined here as the faces of the planar street 
network). The effect of Haussmann modifications on the geometrical 
structure of blocks can be quantitatively measured by the distribution 
of the shape factor <j> (see Methods) shown in Fig. 7. We see that 
before the Haussmann modifications, the distribution of 0 is stable 
and is essentially centered around </> = 0.5 which corresponds to 
rectangles. From 1888, the distribution is however much flatter 




4> 



Figure 7 | Probability distribution of the </> shape factor for the blocks at 
different years. Until 1836, this distribution is stable and we observe a 
dramatical change during the Haussmann period with a larger abundance 
of blocks with small value of (/). These small values correspond to elongated 
rectangle or triangles created by streets crossing the existing geometry at 
various angles. 



SCIENTIFIC REPORTS | 3 : 2153 | DOI: 1 0.1 038/srep021 53 



5 



showing a larger diversity of shapes. In particular, we see that for 
small values of 0 < 0.25 there is an important increase of P{(j>) 
demonstrating an abundance of elongated shapes (triangles and rec- 
tangles mostly) created by Haussmann's works. These effects can be 
confirmed by observing the angle distribution of roads shown on 
Fig. 8 where we represent on a polar plot r(0) = P(0) with P(0) the 
probability that a road segment makes an angle 0 with the horizontal 
line. Before Haussmann's modifications, the distribution has two 
clear peaks corresponding to perpendicular streets and in 1888 we 
indeed observe a more uniform distribution with a large proportion 
of various angles such as diagonals. 

Discussion 

In this paper, we have studied the evolution of the street network of 
the city of Paris. This case is particularly interesting as Paris experi- 
enced large modifications in the 19th century (the Haussmann per- 
iod) allowing us to try to quantity the effect of central planning. Our 
results for central Paris reveal that most indicators follow a smooth 
evolution, dominated by a densification process, despite the import- 
ant perturbation that happened during Haussmann. In our results, 
the important quantitative signature of central planning is the spatial 
reorganization of the most central nodes, in contrast with other 
regions where self- organization dominated and which didn't experi- 
ence such a large-scale structure modification. This structural reor- 
ganization was obtained by the creation at a large scale of new roads 
and avenues (and the destruction of older roads) which do not follow 
the constraints of the existing geometry. These new roads do not 
follow the densification/exploration process but appear at various 
angles and intersect with many other existing roads. 

While the natural, self- organized evolution of roads seems in 
general to be local in space, the Haussmann modifications happen 
during a relatively short time and at a large spatial scale by connect- 
ing important nodes which are far away in the network. Following 
the Haussmann interventions, the natural processes take over on 
the modified substrate. It is unclear at this stage if Haussmann 
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Figure 8 | Radial representation of the angle distribution of road 
segments for 1789, 1826, 1836, 1888. The radial distance r in this plot 
represents the probability to observe a street with angle 6: r = P(6) with 6 
€= [ — 7i/2, nil] and P(6) is the probability to observe an oriented road with 
angle 6 with the horizontal line (see first panel, top left). Until 1836, the 
distribution is peaked around two values separated by approximately 90 
degrees and in 1888, we observe an important fraction of diagonals and 
other lines at intermediate angles. 



modifications were optimal and more importantly, if they were at 
a certain point inevitable and would have happened anyway (due to 
the high level of congestion for example). More work, with more data 
on a larger spatial scale are probably needed to study these important 
questions. 

Methods 

Temporal network data. We denote by G t = G(V t , E t ) the obtained primal graph at 
time t, where V t and E t are respectively the set of nodes and links at time t. The number 
of nodes at time t is then N(t) = \V t \ and the number of links is E(t) = \E t \. Using 
common definitions, we thus have V t = V t - 1 U AV t and E t = E t - 1 U AE P where AV t 
and AE t are respectively the new street junctions and the new streets added in time ] t 
— 1, t] to the network existing at time t — 1. 

The networks for 1789, 1826, 1836, 1888 are extracted from the following maps: 

• 1789: Map of the city of Paris with its new enclosure. Geometrically based on the 
'meridienne de l'Observatoire' and surveyed by Edme Verniquet. Achieved in 
1791. 

• 1826: Road map of Paris surveyed by Charles Picquet, geographer for the King 
and the duke of Orleans. 

• 1836: Cadastre of Paris, Philibert Vasserot. Map constructed according blocks 
and classified according to old districts. 24 Atlas, 1810-1836. 

• 1888: Atlas of the 20 districts of Paris, surveyed by M. Alphand, and L. Fauve, 
under the administration of the prefect E. Poubelle, Paris, 1888. 

All these maps were digitized at the LaDeHiS under the supervision of Maurizio 
Gribaudi, in the framework of a research on the social and architectural transfor- 
mations of parisian neighborhoods between the 18th and 19th centuries. The network 
(and the block structure of figure 1) extracted from the Vasserot cadastre was initiated 
by Anne-Laure Bethe for the program Alpage 29 . 

The networks of 1999 and 2010 are coming from the french Geographical National 
Institute (IGN) on the basis of modern surveys. 

Average route distance. For a network, the shortest path between two nodes is 
defined as the path with the minimum number of links connecting the two nodes. For 
spatial networks, it makes more sense to weight the links with their length: to each 
edge e we thus associate a weight d e defined as the euclidean length of e. We can then 
compute the length I of a path P 

The shortest weighted path is then the one with the minimum total length. The 
average shortest weighted path is also called the average route distance d R . It indicates 
on average how many kilometers you have to walk from one point to the other in this 
spatial network. For a two dimensional network, it is expected 22 that it scales as 

d R ~N 112 (3) 

for a network of size N. In order to compare networks with different numbers of nodes 
AT, it is then natural to compare the rescaled average route distance d R / y/N. 

Betweenness centrality, impact. The nature of the growth process can be 
quantitatively characterised by looking at the centrality of streets. Among the various 
centrality indices available for spatial networks we use here the betweenness centrality 
(BC) 15,17 ' 27 , which is one of the measures of centrality commonly adopted to quantify 
the importance of a node or a link in a graph. Given the graph G t = G(V t , E t ) at time t, 
the BC of a link e is defined as: 



ieV jeV °n 



(4) 



where <7y is the number of shortest paths from node i to node j, while o { fe) is the 
number of such shortest paths which contain the link e. The quantity g(e) essentially 
measures the number of times a link is used in the shortest paths connecting any pair 
of nodes in the network, and is thus a measure of the contribution of a link in the 
organisation of flows in the network. The BC of a node is defined in a similar way 



*(0 



s,teV 



(5) 



where a st {i) denotes here the number of shortest path from node s to t going through 
the node i. 

In order to evaluate the impact of a new link on the overall distribution of the 
betweenness centrality we use the betweenness centrality impact defined in 24 . In the 
graph at time t, we first compute the average betweenness centrality of all the links of 

Gf as: 

(«) 

' eeE t 



(N(t)-l)W)-2) 



where g(e) is the betweenness centrality of the edge e in the graph G t . Then, for each 
link e* £ AE t , i.e. for each newly added link in the time window ]t — 1, t] we consider 
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the new graph obtained by removing the link e* from G t and we denote this graph as 
G t \{e*}. We compute again the average edge betweenness centrality, this time for the 
graph G t \{e*}. Finally, the impact d(e*) of edge e* on the betweenness centrality of the 
network at time t is defined as 



S(e*)- 



E(G,)-g(G,\K»1 



(7) 



The BC impact is thus the relative variation of the graph average betweenness due to 
the removal of the link e*. 

Form factor. The shape or form factor $ of blocks is defined as the ratio of the area of 
the block and the area of the circumscribed circle of diameter D (see 16,23 ) 

A 



kD 2 /4 

The more anisotropic the block and the smaller the factor (j). 



(8) 
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